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ABSTRACT 


Physical  parameters  and  analytical  performance  are  determined  for  an 
analytical  ICP  operated  at  148  MHz,  a  frequency  nearly  three  times  higher 
than  any  previously  reported.  The  iron( I  )  excitation  temperatures  are 

approximately  1.5  times  lower  and  the  electron  densities  are  5  times  lower 


than  at  27 

MHz.  The 

consequences  of  these  changes 

are  lower 

analyte 

and 

background 

continuum 

emission  intensities,  such 

that  the 

signal 

to 

background  ratios  are  decreased  at  the  higher  frequency.  Freedom  from 
interf erences  and  working  curve  linearity  are  unaffected  while  ease  of 
sample  introduction  is  improved.  A  shift  towards  atomic  emission  indicates 
a  deviation  farther  from  LTE  at  148  MHz.  These  effects  are  attributed  to 
the  decrease  in  skin  depth  with  increasing  frequency. 
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1.  INTRODUCTION 

Much  progress  has  been  made  in  the  field  of  ICP-AE5  in  the  two  decades 
since  REED’s  first  reports  [1,21.  Plasmas  have  been  operated  at  a  wide 
range  of  frequencies  since  his  work  at  4  MH2 ,  with  most  systems  today 
operating  at  a  frequency  of  27  MHz .  The  recent  trend  towards  higrer 
frequency  operation  of  ICPs  is  due  to  a  change  in  the  energy  deposition 
region  with  frequency.  It  is  useful  to  think  of  the  energy  deposition 
region  in  terms  of  the  skin  depth,  which  is  defined  as  the  distance  required 
for  the  induced  field  strength  to  fall  off  to  1/e  of  its  value  at  the 
surface,  and  is  inversely  proportional  to  the  square  root  of  the  discharge 
frequency  [31.  All  other  factors  being  equal,  increasing  the  generator 
frequency  moves  the  primary  energy  deposition  region  farther  towards  the 
outside  of  the  plasma.  One  of  the  advantages  of  a  reduced  skin  depth  is 
that  the  analyte,  which  is  injected  along  the  central  channel,  is  less 
likely  to  perturb  the  energy  deposition  process,  and  sample  introduction  is 
therefore  easier  at  higher  frequencies.  The  tradeoff  13  that  the  analyte 
consequently  experiences  a  lower  energy  environment.  CAPELLE  §t_  al .  [4] 
have  investigated  the  influence  of  generator  frequency  at  5,  27,  40,  50  and 
SB  MHz  and  found  that  analytical  performance  is  improved  at  56  MHz  compared 
to  27  MHz.  GUNTER  at  al .  [5]  compared  3,  27  and  50  MHz  ICPs  and  observed 
similar  trends  in  excitation  temperature  and  electron  number  density.  This 
manuscript  charac t er 1 zes  the  operation  of  an  ICP  at  148  MHz,  nearly  three 
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times  the  highest  frequency  investigated  in  previous  studies,  affording  an 
extended  view  of  the  effect  of  generator  frequency  on  the  plasma. 

The  physical  parameters  of  excitation  temperature,  electron  number 
density,  and  ion-to-atom  line  intensity  ratios  have  been  determined  at 
several  power  levels  for  the  two  frequencies  of  27  MHz  and  148  MHz.  A 
number  of  potential  variables  are  held  constant  by  using  the  same  optical 
system  and  thermcmetnc  species  at  each  frequency,  allowing  a  meaningful 
comparison  of  the  two  ICFs.  In  order  to  determine  the  consequences  of  the 
changes  in  the  physical  parameters,  analytical  performance  has  been 
evaluated  for  several  elements  in  terms  of  s ignal-to-background  ratios, 
detection  limits,  calibration  curves,  and  effect  of  interferences.  Ease  of 
organic  sample  introduction  is  also  evaluated. 


2.  APPARATUS  and  TECHNIQUES 


(All  of  section  2  may  be  printed  in  smaller  type) 


2.1.  Instrumentation 

The  27  MHz  RF  generator,  matching  network,  and  standard  size  torch  were 
identical  to  the  ones  previously  described  16,7],  While  the  27  MHz 
generator  allows  a  continuous  adjustment  of  power  level,  the  148. 26  MHz 
crystal-control  led  RF  generator,  a  Hughes  (Fullerton,  CA )  model  HC-300B  UHF 
Command  Transmitter,  provides  a  more  limited  cnoice  of  six  switch-selected 
output  power  levels  of  250.  -100,  550,  900.  12  30,  1500,  ana  2300  watts.  The 
matching  network  was  modified  to  operate  at  i-iS  MH;  py  -epi3cing  the  2—1  2 


turn  load  coil  with  a  one-turn  coil,  and  the  output  pi-network  vacuum 
variable  capacitor  was  replaced  with  a  variable  air-gap  capacitor  machined 
from  S  mm  thick  copper  plates  80  mm  in  diameter,  spaced  by  approximately 
IS  mm . 

The  same  torch  was  used  in  both  systems,  and  was  mounted  in  a  shielded 
enclosure  on  two  Uelmex,  Inc.  (East  Bloomfield,  NY)  dovetail  slides  driven 

by  Rapid-Syn  (Oceanside,  CA  )  stepper  motors  controlled  by  a  Denco  Research, 


Inc.  (Tucson , 

AZ  ) 

SM-2A 

control ler , 

This  allows 

computer-controlled 

positioning  of 

the 

torch 

vertical  1 y 

and  laterally 

with  respect  to  the 

optical  axis. 

A  75 

mm  diameter  x  120 

mm  focal  length 

quartz  lens  was  used 

to  image  the  plasma  with  1:1  magnification  onto  the  entrance  slit  of  a  GCA- 
McPherson  (Acton,  MA  )  EU-700  monochromator  equipped  with  a  Hamamatsu 
(Middlesex,  NJ  >  R212  photomultiplier  tube.  The  electrometer ,  analog-to- 
digital  converter,  and  computer  system  are  similar  to  those  described  by 
HEINE  et.  al .  C71.  All  data  were  taken  under  computer  control  utilizing 
routines  written  in  CONFERS  C81. 

The  nebulizer  is  of  cross-flow  design  constructed  in  our  laboratory, 
with  a  sample  uptake  rate  of  1.3  mL/min  and  an  aerosol  output  of  1.1  L/min. 
A  variable  vent  is  placed  on  the  spray  chamber  to  regulate  the  actual  flow 
to  the  plasma.  A  nebulizer  gas  flow  of  0.7  L/min  was  used  for  both 
frequencies.  Coolant  gas  flow  was  approximately  12  L/mm  and  no  auxiliary 


gas  was  needed. 


Iron  was  chosen  as  the  primary  thermometric  species  for  three  reasons: 

(1)  many  spectral  lines  are  available  in  a  small  wavelength  region,  so 

correction  for  the  variation  of  instrumental  response  with  wavelength  is  not 

necessary,  (2)  the  transition  probabilities  for  iron  are  generally  more 

reliable,  and  (3)  many  previous  studies  of  the  ICP  have  used  iron  so  direct 

comparison  of  results  is  possible.  Excitation  temperatures  were  determined 

by  the  slope  method  [9,101  with  a  linear — 1  eas t -squares  fit,  m  which 

ln(g  A  , / I7i  >  is  plotted  versus  excitation  energy  of  the  spectral  lines,  and 
u  ul 

the  slope  is  equal  to  1/kT.  For  a  transition  from  upper  level  u  to  lower 
level  1,  I  is  spectral  line  intensity,  A  is  the  wavelength,  g^  and  g^  are 

the  statistical  weights  of  each  level,  and  A  ,  is  the  transition 

ul 

probability  for  spontaneous  emission.  Transition  probability  data  are  often 

reported  as  log,_(g,f,  >,  where  f,  is  the  oscillator  strength  and  is 
10  1  lu  lu 

related  to  the  transition  probability  by: 

g,  f,  »  1,499  x  10  ^  A^g  ft  ,  <1> 

1  lu  u  ul 

when  A  j  is  in  sec  *  and  7\  is  in  nm.  Hence  ln(  g^  f  ^/I  >  can  be  plotted 

instead  of  ln<g  A  ,/lT*)  if  log(g,f,  >  values  are  available.  I,  A,  and  g,  f, 
u  ul  1  lu  1  lu 

may  be  on  relative  scales,  but  if  E  is  expressed  in  cm  *  and  T  in  Kelvin, 

exc 

the  slope  will  be  1.4394/T.  The  relative  standard  deviation  of  the  slope  of 
the  line  can  be  used  to  calculate  an  estimated  error,  £T .  The  spectroscop  ic 
constants  of  the  lines  used  are  shown  in  Table  i. 
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Many  of  the  earlier  studies  of  excitation  temperature  in  the  ICP  have 
utilised  the  transition  probability  data  of  CORLISS  and  SQZMAN  [111  or 
CORLISS  and  TECH  C 1 2  3 ,  but  a  systematic  error  in  this  compilation  generally 
leads  to  lower  values  of  excitation  temperature  than  the  corrected  data  C 1 3  3 
More  recent  transition  probabilities  from  BRIDGES  and  KORNBLITH  [141  have 
been  judged  to  be  more  accurate  and  the  data  from  BLACKWELL  et.  al .  [15-211 
to  be  most  accurate  in  a  critical  compilation  of  iron  transition  probability 
data  [221.  Another  commonly  used  set  of  transition  probab  1 1 1 1 les  is  from 
REIF  et.  al .  [231.  The  least  uncertainty  &T  is  obtained  using  the  critically 
evaluated  transition  probabilities  1221  with  the  data  obtained  in  the 
present  study.  Slightly  higher  excitation  temperatures  but  larger 
uncertaint  ies  are  obtained  using  transition  probability  values  from  REIF  et 
al  .  or  BRIDGES  and  KORNBLITH,  while  larger  uncertaint ies  and  markedly  lower 
temperatures  are  obtained  using  the  older  values  [111.  This  is  in  agreement 
with  another  recent  study  of  temperatures  in  the  ICP  using  a  Fourier 
transform  spectrometer  by  FAIRE5  &t  al .  [241. 

Lateral  profiles  of  the  six  wavelengths  taken  at  several  observation 
heights  were  Abel  inverted  C 2 5 3  and  the  inverted  intensities  used  to 
calculate  new  excitation  temperatures .  The  inverted  temperatures  were 
within  the  same  error  bars  in  each  case,  so  Abel  inversion  was  deemed 
unnecessary  for  injected  elements  viewed  on  axis,  in  agreement  with  an 


ear  1 ler  study  C 10  3 


F  i  gure 


shows  a  comparison  of  the  results  obtained  from  a 


representative  data  set  of  the  Fe< I  >  lines  analyzed  with  the  various 
different  sets  of  transition  prooabilities.  The  error  bars  on  the  p i o t 


<  .  «  — %  •  ».*•  .  '•  J  -  _“»•  ,  *"*  »  “  m  *  ••  "  X*  «  *  *  •  ’  *  ’  •  , 

«  4a  *  •  K  a  A  11  'I  <  I  «  .  -  »  *4  a  1  M  1  f  A  T  «  .  »  T  A.  “  •  a*  -f  a  '  »  »  '  -  ’  a*  %'*  ■  “  *  *   ‘  a 


correspond  to  the  uncertainties  in  the  temperatures  calculated  at  each 
observation  height.  These  results  are  very  comparable  in  magnitude  to  those 
presented  in  the  study  by  FAIRES  ejL  al  ■  C  2  4  3  when  using  the  same  transition 
probab i 1  i  t i es  . 

As  a  check  on  the  iron  excitation  temperatures ,  six  lines  of  titanium 
ion  emission  were  also  used.  Spectroscopic  data  for  these  lines  are  shown 
in  Tacle  Z  .  Transition  proaaniiities  are  less  accurately  known  for 
titanium,  so  there  is  mere  uncertainty  in  tnese  resuits.  A  critical 
compilation  C2S]  indicated  better  transition  probabilities  from  the  work  of 
ROBERTS  e_t.  al .  127,281  but  a  much  better  fit  is  obtained  with  the  data  of 
WARNER  1291.  This  again  is  in  agreem_nt  with  an  earlier  study  [101. 

A  similar  excitation  temperature  is  obtained  with  the  titanium  ion  as 
with  the  iron  atom  emission,  as  can  be  seen  in  Fig.  2  for  two  different 
power  levels.  The  higher  titanium  temperatures  at  the  lower  observation 
heights  may  be  due  to  the  much  lower  emission  there  coupled  with  the  much 
higher  sloping  spectral  background,  giving  anomalous  peak  heights  in  the 
data  collection  routine. 
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discharge  even  at  low  concentrations,  and  did  not  enhance  the  H  emission. 

A  silt  width  of  25  microns  was  used  at  27  MHr ,  while  35  microns  was  used  at 

148  MHr  to  compensate  for  the  weaker  emission  at  the  higher  frequency, 

giving  instrumental  half-widths  of  2.05  and  0.07  nm,  respectively.  The 

Doppler  broadening  contribution  can  be  calculated  from  the  relat lonship • 

0  -7  i  /  n 

A>  ,  =  7. 17  :<  10  X  <  T/m  ;  '  2  i 

vre-e  -  ^  =  i  r,  and  m  in  atom;,;  .j  r  1 1  s  [301.  ~  n  i  s  is  -  c  1  e  r  ,•  sensitive  ; 
temperature,  ana  a  half-width  of  0.025  nm  calculated  frcm  a  temperature  of 
5000  K  was  assumed  for  all  measurements.  These  half-widths  when  deconvolved 
from  the  exper i men t a  1 i y  observed  half-width  then  give  the  Stark  half-width, 

which  is  related  to  the  electron  density,  n  ,  by- 
1  /  Z  6 

_  -  2.50  x  10"10  *  n2/3  •:  3  > 

l  /  u  l  /  «i  e 

S  -3 

where  AX^.,  is  in  nm  and  n^  is  in  cm  [301.  The  half-width  parameter  ^  _ 

can  be  obtained  from  tabulations  of  the  theoretical  line  profile 
[31-331.  Since  ^  15  somewhat  dependent  on  temperature  and  electron 

density,  an  iterative  technique  has  most  often  been  employed  110,34,351. 
Uith  the  data  from  the  present  study,  the  error  caused  by  an  order  of 
magnitude  change  in  assumed  electron  density  or  temperature  was  only  4 
percent  compared  to  10  percent  from  the  uncertainty  m  determining  the  half¬ 
widths,  so  the  use  of  a  single  relationship  as  plotted  by  wIESt  1301  is 
justified.  The  relationship  between  observed  n  a i f - width  anc  elec:ron 
density  is  quantified  as ■ 


leg  n  )  =  1.4ol  log,  v  AX  )  +  16.000  (4) 

1 0  e  10  1/2 

when  Ai\  15  expressed  in  nm .  The  overall  uncertainty  in  n  is  estimated 
1/4  e 

to  be  less  than  20  percent.  Based  on  the  results  of  the  excitation 

temperature  study,  Abel  inversions  were  not  deemed  necessary. 

2.4.  6:  oral  to  5  ac/  or  cure  3atios 

Signal  to  oackground  ratios  were  determined  for  several  elements  .res'- 
compromise  conditions  of  1000  W  at  27  MHz  and  1200  W  at  148  MHz.  Two 
analog-to-digital  converters  were  used  to  obtain  the  S/B  ratios  because  of 
the  wide  dynamic  range  involved.  One  ADC  digitized  the  total  analyte-plus- 
bacKground  signal  above  a  zero  set  with  the  light  beam  blocked.  The  other, 
set  to  ten  times  the  sensitivity  of  the  first,  digitized  the  background 
signal  at  a  nearby  wavelength  previously  determined  to  nave  the  same 
intensity  as  at  the  analyte  line.  This  procedure  was  then  repeated  for 
the  range  of  observation  heights. 


2 . S .  Theoretical  LTE  Temperatures 

The  Saha-Eggert  relationship  may  be  used  to  describe  the  equilibrium 

between  the  successive  ionization  stages  of  an  analyte  or  the  argon  support 

gas.  Assuming  Local  Thermal  Equilibrium,  a  unique  n  may  be  calculated  from 

e 

each  temperature  after  the  manner  of  BLADES  l  3S  3 .  Reversing  the  theoretical 

cormulas  to  obtain  the  corresponding  T  .  from  a  measured  n  is  difficult 

n  ,  l  :  e  e 

e 

out  a  thirc-craer  polynomial  was  * ound  to  fit  me  two  parameters  well  z  *  er 


the  range  of  T  3  6500  to  10,000  K  and  n  =  1  x  10  to  1  x  10  cm.  This 

e 

relationship  is: 

T  .  _  -  S475  +  1026 ( X  )  -  64.32(X2)  +  53.12(X3  )  (5) 

n  ,LTE 
e 

where  X  -  log,„(n  >  -  13  in  order  to  reduce  round-off  errors.  This  formula 
10  e 

was  useful  in  determining  theoretical  equilibrium  temperatures  based  on 
experimental  electron  density  measurements. 

3.  RESULTS  and  DISCUSSION 

3.1.  Physical  Parameters  of  the  Oischarne 

3.1.1.  Excitation  Temperatures .  Figures  3  and  4  present  iron 
excitation  temperatures  determined  at  the  several  power  levels  as  a  function 
of  observation  height  in  the  plasma,  using  the  critically  evaluated 
transition  probabilities.  When  the  different  transition  probability  data 
sets  and  power  levels  are  taken  into  account,  a  real  comparison  can  be  made 
between  this  study  and  those  of  other  researchers  who  have  investigated  the 
vertical  spatial  profile  of  the  excitation  temperature  at  27  MHz.  FURUTA 
and  HORLICK  [37]  made  observations  at  a  single  fairly  high  power  level  of 
1800  W,  which  would  be  expected  to  give  high  temperatures ,  but  used  CORLISS 
and  BOZMAN’s  till  transition  probabilities,  which  then  bring  the  reported 
temperatures  down.  KAWAGUCHI  j*_t  al  .  [38]  made  observations  at  three 
different  power  levels  and  usea  transition  probab 1 1  i t i es  from  REIF  e_t.  al  ■ 
[231,  which  as  expected  give  slightly  higner  temperatures  than  obtained  m 


this  study . 


Both  of  these  groups  observed  a  similar  dip  in  excitation 
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temperature  at  low  observation  heights,  and  this  dip  is  seen  to  move  lower 
in  the  plasma  as  power  is  increased  in  the  present  investigation  as  well. 

While  the  27  MHz  plasma  exhibits  a  strong  power  level  dependence  of  the 
excitation  temperatures ,  the  very  high  frequency  plasma  does  not.  This  is 
in  accordance  with  their  visual  appearance  when  a  1000  ppm  solution  of 
yttrium  is  aspirated.  With  a  nebulizer  gas  flow  of  0.7  L/min,  the  blue 
ionic  emission  region  m  the  so-called  "Normal  Analytical  Zone' 
[39, *10]  begins  approximately  15  mm  above  the  load  coil  at  the  lower 
frequency,  with  increasing  intensity  observed  with  increasing  power.  The 
VHF  discharge  at  900  W  does  not  develop  this  emission  region,  but  remains  a 
red  color  throughout  in  the  standard  size  torch.  At  higher  powers,  a  small 
blue  patch  develops  about  15  mm  above  the  load  coil,  but  remains  very  faint 
even  at  1500  W. 

That  higher  frequencies  yield  lower  excitation  temperatures  in  the  ICF 
has  been  previously  demonstrated  [4,53,  so  a  lower  temperature  is  not 
unexpected  for  the  148  MHz  plasma.  The  effects  of  different  power  levels 
and  specific  observation  heights  were  not  addressed  in  their  paper,  so  the 
linear  relationship  between  temperature  and  frequency  shown  by  CAPELLE  et 
al .  [4]  may  only  be  fortuitous.  If  that  relationship  is  extended  to 
148  MHz,  negative  temperatures  would  have  been  expected,  clearly  an 
unrealistic  sitation.  Since  calcium  emission  lines  were  used  for  the 
excitation  temperatures  in  the  other  frequency  study  [53,  only  the 
similarity  in  trend  can  be  noted  there. 


3.1.2.  Electron  Densities.  The  calculated  electron  densities  as  a 
function  of  power  and  observation  height  in  the  two  plasmas  are  shown  in 
Figs.  5  and  6.  Many  researchers  have  reported  electron  densities  only  at  a 
single  observation  height  and  power  level.  For  example,  the  value  of 
15  -3 

5  x  10  cm  at  1500  W  and  15  mm  observation  height  determined  by  UCHIDA  e_t. 

15  -3 

al .  [41]  at  27  MHz  is  slightly  higher  than  our  value  of  3.5  x  10  cm 

under  similar  conditions,  wmle  1.6  x  10^  cm  0  from  REZAAIYAAN  and  MIEFTJE 
[35]  is  slightly  lower.  ALDER  et.  a  1 .  [42]  showed  electron  densities  over 
the  range  10  mm  -  25  mm  observation  height  for  one  power  of  1200  W  which  are 
only  slightly  higher  than  the  observed  curve  for  1250  W  at  27  MHz.  Axial 
values  extracted  from  the  spatial  data  of  CAUGLIN  and  BLADES  [45]  are  also 
similar.  Although  these  spatial  data  indicate  higher  electron  densities 
outside  of  the  central  channel  at  low  observation  heights,  the  lack  of  Abel 
inversion  techniques  should  not  affect  the  validity  of  the  lateral 
measurements,  especially  for  heights  greater  than  8  mm,  because  the  electron 
densities  in  the  present  study  were  derived  from  H|t,  emission  from  injected 

water,  which  is  mainly  confined  to  the  anlayte  channel. 

Figures  5  and  5  show  a  fairly  smooth  drop  in  electron  density  with 
observation  height  and  a  direct  dependence  on  power  level  for  both  ICPs. 
The  148  MHz  discharge  exhibits  about  half  an  order  of  magnitude  lower 
electron  density  overall,  and  in  contrast  to  the  excitation  temperature,  a 
direct  power  level  dependence  is  observed.  The  lower  electron  density  is 

expected  from  other  workers  ’  results  at  higher  frequencies  [10,44]  but  again 
linear  trend  shown  by  CAFELLE  s_t_  a  1 .  [4]  15  continued  toe  preoicteo 


if  the 


electron  density  at  14-8  MHz  would  be  seven  orders  of  magnitude  lower  than  is 
actually  observed. 


3.1.3.  Ion  to  Atom  Line  Ratios.  It  has  been  stated  that  a  relative 
idea  of  electron  density  can  be  easily  obtained  from  the  ratio  of  ionic  to 
atomic  emission  from  a  suitable  element  C41.  Magnesium  is  often  used 

because  of  the  proximity  of  the  Mgi  I  )  line  at  235.  Zi  nm  and  two  Mgi II'  iines 

at  273.55  and  230.27  nm .  When  l_7E  is  assumed,  7  =7  =  7  __  and  the 

e  :<  c  ion  Lit 

Saha  and  Boltzmann  equations  can  be  combined  to  yield  the  following 
relationship  between  electron  density  and  ion/atom  emission  ratios; 


T-v  2ir  < 2  m  k  )3/^  g+A-l’/-v-» 

I_  . _ S _  _ _ L 

1°  h3  g’^/y 


exp 


-E  -  E  + 
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kT 


LTE 
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where  n  =  electron  density  ( in  m  ), 
e 

-31 

m  =  electron  mass  <9.11  x  10  kg), 
e 

-23  -1 

k  *  Boltzmann’s  constant  (1.381  x  10  J/K  or  0.6947  cm  /K), 

-34 

h  *  Planck's  constant  <6.626  x  10  J-s), 

I°,I+  »  emission  intensity  of  neutral  and  ionized  species, 

g’,g+  =  statistical  weight  of  the  emitting  levei  in  the  atom  ana  ion, 

4°,nJ'  =  transition  probability, 


V , V  =  wavelength  of  emission. 


E  *  ionization  energy, 
ion 


r' 

Er  ,E  =  excitation  energy  of  the  emitting  level 
ex  c  exc 


in  the  atom  and  ion, 

*  temperature. 

The  emission  intensity  ratio  is  dependent  not  only  on  electron  density,  but 

also  or  the  temperature.  Since  the  excitation  temperature  has  been  shown  :: 

vary  dramatically  wit n  power  and  observation  height,  a  simple  measurement  or 

the  Mgt 1 1  )/Mg( I >  ratio  will  not  strictly  follow  the  electron  density,  as  can 

be  seen  in  Figs.  7  and  8,  which  show  the  observed  Mg  ratios  for  the  two 

frequencies.  Electron  density  values  calculated  from  observed  lon/atom 

intensity  ratios  and  excitation  temperatures  are  often  orders  of  magnitude 

lower  than  n  values  obtained  from  Stark  broadening  measurements  [  34  3. 
c 

Indeed,  experimental  lon/atom  emission  intensity  ratios  have  been  proposed 
as  a  measure  of  the  deviation  from  LTE  1433. 

The  values  reported  for  the  ratio  of  Mg(II)/Mg(I)  emission  at  single 
observation  height?  *nd  power  levels  by  KALNICKY  et.  al .  1343  and  CAPELLE 
et  al .  [43  are  about  twice  the  observed  value  of  thi3  ratio  under  similar 

conditions,  but  the  vertical  spatial  behavior  shown  by  FURUTA  and  HORLICK 
[373  agrees  well  with  these  observations.  The  vertical  spatial  profile 
deduced  from  data  presented  by  CAUSHLIN  3nd  BLADES  [433  is  also  similar  in 
shape  to  that  obtained  in  the  present  study.  The  differences  may  be  due  to 
the  lack  of  Abel  inversion  of  the  ion  and  atom  intensities  here. 

A  dramatic  effect  of  observation  height  is  seen  at  37  MHz,  roughly 


zones  of  the  plasma  observed  with  yttrium 


corresponding 


to  the  different 


k 


k 


t- 
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introduced:  lower  in  the  plasma  principally  atomic  emission  is  observed, 
while  the  ionic  emission  is  maximized  in  the  Normal  Analytical  Zone  and  then 
falls  off  with  increasing  height.  In  the  148  MHz  plasma,  the  blue  ionic 
emission  is  much  less  intense,  and  the  Mg(II>/Mg<I)  ratios  are  lower  as 
well.  In  addition,  the  influence  of  power  level  is  less  pronounced.  These 
effects  are  consistent  with  the  changes  in  excitation  temperature  and 
electron  censity  seen  aoove. 

3.2.  Analytical  Performance 

3.2.1.  Signal  to  Background  Ratios.  Generally,  the  signal  to 
background  ratio  was  found  to  be  strongly  dependent  on  observation  height. 
Typical  behavior  is  shown  in  Fig.  9  for  copper( I  )  at  324.75  nm.  Low  in  the 
plasma,  the  S/B  ratio  was  low,  and  increased  to  a  maximum  about  30-33  mm 
above  the  load  coil.  Applying  the  method  of  WINGE  e£  al ■  C45J,  one  might 
conclude  that  this  observation  height  would  provide  the  best  detection 
limit.  The  best  actual  detection  limits  however  are  not  strictly  a  function 
of  the  S/B  ratio  [3,463.  At  this  height,  the  background  was  indeed  very 
low,  but  so  was  the  net  signal,  and  with  a  limited  amount  of  gain  available 
in  the  electrometer ,  the  lowest  detection  limits  were  not  obtained.  The 
maximum  net  signal  for  most  elements  was  found  to  be  17-20  mm  above  the  load 
coil  for  both  frequencies.  In  this  spatial  region,  the  actual  detection 
limits  were  reasonably  close  to  those  approximated  from  the  S/B  ratios. 
Results  for  the  observation  heights  at  which  the  detection  limits  were 
obtained  are  summarized  and  compared  in  Tadie  3.  The  signal  to  background 
ratio  decreases  with  increasing  frequency  in  ail  cases  except  Ca<  I  at 
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4Z2.67  nm.  Because  the  background  intensity  decreases  similarly  for  ail  the 
transitions,  the  3/B  ratio  changes  primarily  in  response  to  the  decreasing 
net  signal  intensity,  with  ion  lines  more  sensitive  to  the  frequency  change 
than  the  atom  lines. 

3.2.2.  Detection  Limits.  A  compromise  observation  height  of  19.5  mm 
above  the  ioad  coil  based  on  the  signal  to  background  ratio  measurements  was 
used  for  ail  elements  except  calcium.  The  maximum  net  signal  'or  calcium 
was  found  to  be  13  mm  above  the  load  coil  for  both  frequencies  and, 
interestingly,  for  both  the  atom  and  ion  emission  lines. 

The  detection  limit  was  based  on  the  concentration  of  analyte  that  gave 
an  average  signal  equal  to  twice  the  peak-to-peak  noise  of  the  baseline 
actually  evaluated  at  that  concentration.  A  time  constant  of  0.1  sec  was 
used.  The  results  at  27  MHz  are  generally  5-10  times  higner  than  those 
tabulated  by  WINGE  et.  al .  [45],  but  are  not  unreasonable  considering  the 
differences  in  instrumental  ion . 

Since  identical  procedures  and  instrumentation  were  used  at  both 
frequencies,  a  valid  comparison  can  be  made  between  these  two  plasmas. 
Results  are  summarized  in  Table  4.  In  all  cases  except  the  Fe< I >  and  Ceil) 
lines,  the  detection  limit  at  148  MHz  is  higher  than  that  at  27  MHz.  The 
trend  in  general  follows  the  trend  in  signal  to  background  ratios;  the  ion 
emission  lines  decrease  more  than  the  atom  lines,  so  the  atom  line  detection 
limits  are  closer  to  the  27  MHz  values  than  those  determined  from  the  ion 
lines.  This  is  in  contrast  to  the  results  found  previously  [ 4  ] ,  in  which 
3/B  ratios  ana  detection  limits  improve  with  frequency  up  to  56  MHz. 


Apparently,  analytical  performance  improves  at  twice  the  conventional 
frequency,  but  does  not  continue  the  trend  with  a  six-fold  increase. 

3.2.3.  Calibration  Curves.  Calibration  curves  are  linear  and  extend 
over  four  or  more  orders  of  magnitude  for  both  frequencies,  as  can  be  seen 
from  the  representat i ve  curves  shown  in  Fig.  10.  Calcium  was  the  only 

element  fauna  to  fall  off  significantly  at  1000  ppm  at  Z7  MHz,  out  w as  founa 
to  be  linear  to  1300  ppm  at  148  MHz.  All  other  curves  were  founa  to  be 
strictly  linear  at  both  frequencies,  with  similar  dynamic  ranges. 

3.2.4.  Interf erences .  The  vaporization  interference  of  phosphate  on  a 
5  ^g/mL  solution  of  calcium  was  investigated  at  Doth  frequencies,  and  was 
not  found  to  be  of  significant  magnitude  up  to  molar  ratios  of  700:1 

PO^  -Ca.  At  a  ratio  of  1300:1,  a  5  percent  depression  of  the  atomic 
4 

emission  occurs  for  both  frequencies,  as  shown  in  Fig.  11.  The  ionic 
emission  is  essentially  unaffected  at  27  MHz,  while  at  i48  MHz  a  slightly 
greater  depression  is  seen.  The  six-fold  increase  in  frequency  does  not 
reduce  the  ICP’s  insensitivity  to  this  inter f erence . 

3.2.5.  Qroanics .  In  order  to  determine  the  relative  ease  of  sample 
introduction,  many  common  organic  solvents  were  aspirated  into  the  two 
plasmas,  including  acetone,  benzene,  hexane,  methanol,  and  xylene.  "Organic 
iron"  solutions  were  prepared  from  ferrocene  aissolvea  m  benzene  and 
<yiene.  In  order  to  operate  the  Z7  MHz  ICP  with  the  organic  solvents, 
drastic  retuning  (several  turns  of  eacn  variable  capacitor  in  the  matormg 


network)  and  a  minimum  power  level  of  1500  W  was  necessary.  Less  drastic 
retuning  was  also  necessary  for  different  classes  of  compounds.  Even  under 
these  conditions,  it  was  difficult  to  keep  the  plasma  lit.  At  143  MHz,  only 
the  more  volatile  species  required  a  slight  retuning  (1/8  turn  of  one  of  the 
capacitors)  and  power  levels  as  low  as  650  W  operated  satisfactorily. 

The  difference  in  the  appearance  of  the  two  plasmas  was  striking  when 
benzene  was  aspirated.  At  27  MHz.  the  usual  'green  tongue"  was  seen  m  a 
generally  green  plasma,  with  a  slight  amount  of  yellow  emission  fron 
incomplete  combustion  accompanied  by  some  soot  formation.  At  148  MHz,  the 
yellow  emission  was  more  intense,  giving  an  extremely  bright  plasma  and 
larger  quantities  of  soot.  All  other  solvents  in  both  plasmas  were  simply 
green.  Because  of  this  high  background  emission,  the  most  sensitive  scales 
of  the  electrometer  could  not  be  used  at  the  higher  frequency  with  benzene 
as  the  solvent,  and  xylene  was  used  instead. 

At  27  MHz  and  1500  W,  the  Fe< I  )  line  at  371.99  nm  gave  16  times  more 
signal  than  the  ion  line  at  259.94  nm,  in  contrast  to  the  behavior  with 
aqueous  solutions.  The  atom  line  at  148  MHz  and  1200  LI  provided  the  same 
3ize  signal  as  at  27  MHz,  with  approximately  equal  noise  on  the  same 
electrometer  setting.  Sensitivity  is  therefore  about  equal  for  the  two 
frequencies  with  organic  solvents,  but  the  high  frequency  plasma’s  improved 
tolerance  for  organic  solvents  at  lower  powers  is  a  distinct  advantage. 


It  has  recently  been  suggested  that  the  electron  density,  cemg  a 


single-valued  parameter  amenable  to  determinat ion  without  already  assuming 


Local  Thermal  Equilibrium,  is  most  appropriate  to  characterise  a  discharge 

in  LIE,  rather  than  an  excitation  temperature,  which  may  be  dependent  on  the 

element  involved  and  the  excitation  energies  of  the  levels  used 

[42,43,481.  This  provides  a  basis  for  comparing  the  deviation  from  LTE  in 

an  actual  plasma.  The  observed  n  data  translate  to  the  T  ,  __  curves 

e  n  ,LTE 

e 

shown  in  Figs.  12  and  13.  Transition  probabilities  taken  from  SMITH  and 

LISZT  [49]  and  ether  spectroscopic  data  used  in  the  calculation  of 

theoretical  intensity  ratios  0y  Eqn.  (6)  are  listed  in  Table  5.  The 

parameter  b  introduced  by  CAUGHLIN  and  BLADES  is  then  simply  obtained: 
r 


<1*71°  > 


expt 


(  7  ) 


<  I  V I  *  > 


LTE 


Values  of  b  calculated  from  these  data  are  shown  in  Figs.  14  and  15. 
r 

Since  electron  densities  were  found  to  fail  smoothly  with  height,  tne  b 


values  follow  the  spatial  behavior  of  the  Mg< 1 1  )/Mg(  I  >  intensity  ratios.  In 

all  ■  cases,  b  is  less  than  one,  reflecting  the  use  of  T  .  __  rather  than 
r  n  ,LTE 

e 

T  in  the  LTE  calculation.  The  values  at  27  MHz  are  slightly  less  than 

exc 

those  shown  by  CAUGHLIN  and  BLADES,  largely  due  to  the  difference  in  Mg 
ratios.  It  is  not  clear  at  this  time  what  range  of  variation  may  be 
expected  from  laboratory  to  laboratory  due  to  differences  in  instrumentation 
or  techniques.  Comparing  the  two  frequencies  in  a  single  laboratory,  it  can 
be  seen  that  the  148  MHz  ICP  deviates  muen  more  strongly  from  LTE  tnan  tne 
2"’  MHz  ICP.  As  observation  height  is 
similar  state  of  non-LTE. 


increased,  the  two  plasmas  approach  a 


3.4-.  Pi  5cuss  ion 


The  detection 

limits 

i  n 

the  148 

MHz 

plasma  are  higher  than  those 

a  t 

27  MHz  primarily 

because 

the 

s  igna  1 

has 

fallen  off 

faster  than 

the 

background.  The 

increase 

in 

detection  limits  roughly 

corresponds  to 

t  he 

decrease  in  signal  to  background  ratio.  At  all  wavelengths  observed,  the 
background  was  5-12  times  lower  at  148  MHz.  This  is  primarily  due  to  the  3- 

foid  drop  m  electron  density  from  27  MHz  and  1000  'J  to  14S  MHz  and  1Z30W, 
since  the  background  continuum  emission  is  proportional  to  the  square  of  tne 

electron  density  C441. 

The  decrease  in  electron  density  does  not  scale  directly  with 
frequency,  which  is  understandable  when  the  skin  depth  is  considered.  The 
skin  depth  decreases  by  a  factor  of  2.4  over  this  frequency  range  if  the 
same  power  levels  are  maintained,  and  this  directly  affects  the  relationship 
between  the  energy  deposition  region  and  the  analyte  channel.  Comparing 
equivalent  power  levels  at  each  frequency,  i.e.  1200  W  and  1500  U,  the 
electron  density  is  5  times  lower  at  148  MHz,  which  is  twice  the  change  m 
skin  depth.  It  is  difficult  to  apply  this  hypothesis  to  studies  by  other 
researchers  at  intermediate  frequencies  because  of  the  different 
experimental  conditions  employed.  GUNTER  §i  a  1  .  C  5  3  did  not  directly 
determine  the  electron  densities,  and  CAPELLE  si  a  1 .  C4]  only  show  measured 
electron  densities  for  5,  27  and  40  MHz.  Since  the  5  MHz  ICP  was  operated 
at  five  times  the  power  of  the  other  two,  direct  comparison  cannot  ce  "ace. 
From  27  to  40  MHz,  the  skin  depth  is  expected  to  decrease  by  3  factor  c  *' 
1.2,  wmle  their  oDserved  n^  drops  by  1.7.  Differences  i  r  spp.iea  r  ^ 

other  experimental  parameters  not  reported  m  their  pace'-  mev  ce  r-e  - 


this  result  does  not  compare  directly  with  the  present  study.  Scaling  with 


skin  depth  rather  than  directly  with  frequency  better  explains  the  observed 
magnitude  of  the  decrease  in  electron  density  with  frequency. 

The  reason  behind  the  decrease  in  emission  intensity  is  not  so  clear. 
The  excitation  temperature  drops  by  a  factor  of  1.5  over  this  frequency 
range.  Neither  the  excitation  temperature  nor  the  decrease  in  analyte 
emission  seem  to  be  simply  correiateq  witn  the  orange  in  s r.  i  n  cepth.  'he 
relative  loss  of  net  signal  intensity  might  ae  roughly  -elated  to  t  re 
excitation  energy  of  the  emitting  levels  of  the  atom  lines,  but  this 
correlation  does  not  hold  for  the  ion  lines,  even  when  the  ionization 
potential  is  taken  into  account.  FURUTA  and  HORLICK  [37]  have  proposed  that 
neutral  atom  excitation  is  a  resuit  of  electron  collisions,  while  -oilisions 
with  excited  argon  species  are  responsible  for  tne  generation  and  excitation 
of  ionic  species.  The  lower  electron  densities  would  lead  to  the  lower 
emission  intensities  under  this  scenario,  and  the  different  behavior  of  the 
ion  lines  could  be  rat lona 1 1  zed . 


'Jalues  of  b^_  less  than  one  are  said  to  be  indicative  of  an  ionizing 


non-equilibrium  plasma  undergoing  rapid  heating,  charac t er i zed  by  a  relative 
over-population  of  ground  state  levels  compared  to  LTE  [44],  The  shift 
farther  from  LTE  is  evidenced  by  the  shift  of  the  Mg  lon/atom  intensity 
ratio  towards  atom  emission  with  increased  frequency.  The  same  oenavicr  is 
seen  for  the  other  elements  for  whicn  both  atomic  and  ionic  emission  data 
were  ontamed-  at  I?  MHz  the  Feilli  line  is  approx  imat  ei  y  six  times  the 
intensity  of  the  Fail'  line,  while  it  is  only  aDcut  ha i f  as  intense  at 
1 4  f  v'h; .  Simi  lariv,  at  I~  MHz  the  calcium  ion  line  is  103  times  mere 
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intense  than  the  atom  line,  but  at  14-8  MHz  they  are  of  approximately  equal 
intensity.  Two  effects  have  combined  here  to  cause  this:  the  drop  in 
electron  density  serves  to  shift  ionization  equilibra  away  from  the  ion,  and 
the  greater  relative  over — population  of  the  ground  state  causes  the  atom 
lines  to  decrease  less  relative  to  the  ion  lines  as  frequency  increases. 

4.  CONCLUSION 

These  investigations  of  an  analytical  ICP  at  a  frequency  nearly  three 
times  higher  than  any  previously  reported  demonstrate  the  effects  of  the 
reduced  skin  depth  at  higher  generator  frequencies.  This  produces  a  less 
energetic  environment  character  iced  by  lower  excitation  temperatures  and 
electron  densities,  resulting  in  both  lower  background  and  analyte  emission 
from  the  UHF  plasma.  The  background  intensity  is  proportional  to  the  square 
of  the  electron  density,  which  scales  with  the  skin  depth  rather  than 
directly  with  the  frequency.  The  analyte  emission  is  a  more  complex 
function  which  falls  off  faster  than  the  background  over  this  extended 
frequency  range.  The  observed  shift  away  from  ionic  emission  is  indicative 
of  a  shift  farther  from  LTE  as  skin  depth  is  reduced.  The  ICP’s  immunity  to 
interelement  inter f erences  is  not  lost  at  148  MHz  and  ease  of  sample 
introduction,  especially  of  organic  species,  is  improved. 
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FIGURE  CAPTIONS 


Fig.  1.  Iron(I)  excitation  temperatures  determined  using  several  different 
transition  probability  data  sets  at  27  MHz  and  1000  W.  Transition 
probabilities  from:  |  FUHR  et  [22],  O  BRIDGES  and  KORNBLITH  [14], 
▲  REIF  et  [23],  Q  CORLISS  and  BOZMAN  [11]. 

Fig.  2.  Titanium(II)  and  iron(I)  excitation  temperatures  determined  at 
27  MHz.  |  titanium  at  1000  W,  O  iron  at  1000  W,  A  Titanium  at  1250  W, 
I  1  iron  at  1250  W. 

Fig.  3.  Iron(I)  excitation  temperatures  at  27  MHz.  |  800  W,  O  1000  w* 
A  1250  W,  Q  1500  W. 

Fig.  4.  Iron(I)  excitation  temperatures  at  148  MHz.  Q  900  W,  A  1200  W, 
□  1500  W. 

Fig.  5.  Electron  densities  determined  from  Stark  broadening  at  27  MHz. 

■  300  W,  O  1°00  A  1250  W,  □  1500  W. 

Fig.  6.  Electron  densities  determined  from  Stark  broadening  at  148  MHz. 

O  900  W ,  A  1200  W,  Q  1500  W. 


Fig.  7.  Magnesium  ion  (279.55  ran)  to  atom  (285.21  nm)  emission  intensity 
ratios  determined  at  27  MHz.  |  800  W,  Q  1000  W,  A  1250  W,  Q  1500  W. 

Fig.  8.  Magnesium  ion  (279.55  nm)  to  atom  (285.21  nm)  emission  intensity 
ratios  determined  at  148  MHz.  Q  900  W,  A  1200  W,  Q  1500  W. 

Fig.  9.  Copper(I)  signal  to  background  ratios  and  relative  signal 
intensities  determined  at  324.75  ran  under  compromise  operating  conditions, 
o  27  MHz  and  1000  W,  □  148  MHz  and  1200  W. 

Fig.  10.  Analytical  calibration  curves  for  Q  Fe(II)  (259.94  nm)  and 

O  Cu( I )  (324.75  nm)  at  27  MHz  and  1000  W,  |  Fe(I)  (371.99  nm)  and 

A  Cu(I)  (324.75  ran)  at  148  MHz  and  1200  W. 

Fig.  11.  Effect  of  phosphate  on  calcium  ion  and  atom  emission  intensity 

from  a  3  ug/ml.  solution.  Q  27  MHz  and  1000  W,  Q  148  MHz  and  1200  W. 


Fig.  12.  LTE  temperatures  determined  from  electron  densities  at  27  MHz. 
||  800  W,  O  1000  M»  A  1250  w,  0  1500  W. 

Fig.  13.  LTE  temperatures  determined  from  electron  densities  at  148  MHz. 
Q  900  W,  ▲  1200  W,  Q  1500  W. 

Fig.  14.  Values  of  br  *  (I+/I°)expt  7  calculated  for  27  MHz* 

■  800  W,  O  1000  w»  A  1250  Wf  0  1500  W. 

Fig.  15.  Values  of  bp  *  (1+/I°)expt  /  (i+/iC\te  calcu1ated  for  X48  MHz* 
O  900  W,  A  1200  W,  0  1500  W. 
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